Spin dynamics in the pseudo-gap state of a high-temperature superconductor 
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The pseudogap is one of the most pervasive phenomena of high temperature superconductors 
It is attributed either to incoherent Cooper pairing setting in above the superconducting transition 
temperature T c , or to a hidden order parameter competing with superconductivity. Here we use 
inelastic neutron scattering from underdoped YBa2Cu30e.6 to show that the dispersion relations of 
spin excitations in the superconducting and pseudogap states are qualitatively different. Specifically, 
the extensively studied "hour glass" shape of the magnetic dispersions in the superconducting state 
@, H> El is no longer discernible in the pseudogap state and we observe an unusual "vertical" disper- 
sion with pronounced in-plane anisotropy. The differences between superconducting and pseudogap 
states are thus more profound than generally believed, suggesting a competition between these two 
states. Whereas the high-energy excitations are common to both states and obey the symmetry of 
the copper oxide square lattice, the low-energy excitations in the pseudogap state may be indicative 
of collective fluctuations towards a state with broken orientational symmetry predicted in theoretical 
work [110, !• 

Original reference: Nature Physics 3, 780 (2007). 



The pseudogap (PG) manifests itself in various exper- 
imental probes as a depiction of spectral weight upon 
cooling below a doping-dependent characteristic temper- 
ature T* . The origin of the pseudogap remains the focus 
of significant debate. For instance, recent experiments 
have provided evidence of vortex-like excitations below 
T* akin to those in the superconducting (SC) state 
Other experiments [lfj suggest the presence of a novel 
magnetic order breaking time-reversal symmetry, again 
setting in around T*. 

Inelastic neutron scattering directly probes the 
microscopic magnetic dynamics and can thus serve 
as particularly incisive tests of microscopic models of 
the cuprates. Recent work on non-superconducting 
La2- x (Sr, Ba) x Cu04 and in the SC state of 
YBa2Cu306+x has uncovered tantalizing evidence 
of a "universal" spin excitation spectrum independent of 
material-specific details @, H, 0, U > 12] • The dispersion 
surface comprises upward- and downward-dispersing 
branches merging at the wave vector Qaf, which 
characterizes antiferromagnetic order in the undoped 
parent compounds. The spectrum thus resembles an 
"hour glass" in energy-momentum space. Investigations 
of the SC state of twin-free YBa2Cu306+ x (i.e., single 
crystals with a unique orientation of the two in-plane 
axes a and b throughout the entire volume) have revealed 
a two-dimensional geometry of the lower branch, with 
a modest in-plane anisotropy of the spectral weight 
that increases with decreasing excitation energy [1 31 ] - 
The "hour glass" spectrum and its in-plane anisotropy 
have stimulated theoretical work in the framework of 
both itinerant and local-moment pictures of the electron 
system 0, EE EE EE EE El • 



Previous measurements of twinned YBa2Cu306+i- 
samples appeared to indicate that the spin excitation 
spectrum in the PG state is simply a broadened version 
of the spectrum in the SC state [2^, 21, However, 



the detailed investigation of the geometry and dispersion 
of spin excitations in the PG state reported here reveals 
that both spectra are distinctly different. 

The measurements were performed on an array of twin- 
free single crystals of underdoped YBa 2 Cu306.6 with 
superconducting T c = 61 K (see Methods). Fig. 1 
shows representative raw data in the SC state at 5 K, 
in the PG state just above T c , and at room temperature. 
The incommensurability S and the spectral weight of the 
constant-energy cuts are generally different in the two in- 
plane directions in reciprocal space, H and K, and show 
a complex dependence upon energy and temperature (see 
Methods for definitions). A synopsis of the entire data 
set is shown in Fig. 2 at two temperatures above and 
below T c . 

The central result of this paper is the change of the 
topology of the dispersion surface from the SC to the PG 
state. First, we focus on the spectrum in the SC state 
(Fig. 1, 5 K). Starting from low excitation energies, 5 first 
decreases with increasing energy, so that the incommen- 
surate peaks merge at Qaf at an energy of E Tes = 37.5 
meV. For E > E Ies , S increases again, so that the spec- 
trum forms the "hour-glass" dispersion (Fig. 2a, b) al- 
ready familiar from previous work 

BflEHnii- In the 

PG state, however, the singularity at E Tes is no longer 
discernible, and 5 is only weakly energy-dependent over 
a wide energy range including E ies (Fig. 2c, d and Sup- 
plementary Information). This constitutes a qualitative 
difference of the spectra in the SC and PG states. 
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FIG. 1: Energy evolution of the in-plane magnetic excitations 
around Qaf for different temperatures. In a,b, the energy 
transfer was fixed to 60 meV, in c,d to 47 meV, in e,f to 
37.5 meV and in g,h to 33.5 meV. Panels a,c,e,g show scans 
along the a-axis and panels b,d,f,h scans along the 6-axis. 
The lines are the results of fits to Gaussian profiles. We show 
the raw triple-axis data; the only data processing applied is 
a subtraction of a constant at 250 K and 290 K in order 
to account for the increased background from multi-phonon 
scattering. Corrections for the Bose factor are small and were 
not applied to the data. The final wave vector was fixed to 
2.66 A below 38 meV and to 4.5 A above. Error-bars 
indicate the statistical error. 



The spin excitations in both states also differ markedly 
with respect to their a-6-anisotropy. In the SC state, 
constant-energy cuts of the magnetic spectral weight at 
energies below E res form ellipses with aspect ratios of 
about 0.8 (Figs. lg,h and 3a). In the PG state, the Q- 
extent of the signal decreases significantly along b* such 
that incommensurate peaks can hardly be resolved, Fig. 
lh. In contrast, the flat-top profiles along a* are well 
described by two broad peaks displaced from Qaf , Fig. 
lg, and we can set an upper bound of 0.6 on the ratio of 
5 along b* and a* . The intensity distribution of the high- 
energy excitations, on the other hand, hardly changes 
between the SC and PG states (Figs. la,b and 3b). It 
is also much more isotropic than the low-energy profiles. 
Specifically, the ratio of S along a* and b* is 0.95 ±0.07 as 
obtained from an analysis of the time-of-flight and triple 
axis data. Moreover, there is no amplitude suppression 
along b* and the anisotropy in the peak widths is also 
less pronounced than at low energies (0.16 ± 0.02 and 



0.13 ± 0.02 r.l.u. along a* and b* , respectively). 

Upon cooling below T c , the spectral rearrangement as- 
sociated with the formation of the downward-dispersing 
branch of the "hour-glass" results in a sharp upturn of 
the intensity at points along this branch (Fig. 4a, b), 
while at Qaf and 30 meV there is only a broad max- 
imum at T c (Fig. 4c). This is a further manifestation 
of the qualitative difference between the SC and PG 
state spectra. In the PG state, the spectral weight at 
energies at and below E Tes declines uniformly with in- 
creasing temperature at all Q-values and vanishes around 
T* w 200 K. T* is comparable to the temperature be- 
low which the PG manifests itself in other experimental 
probes 0, 0, [l(|. Corresponding constant-energy cuts 
show that for T > T*, the low-energy spectral weight is 
severely depleted over the entire Brillouin zone below an 
energy E* ~ 40 meV, Fig. le-h (see also Supplemen- 




FIG. 2: Colour representation of the magnetic intensity, ob- 
tained from triple-axis scans. Panels a,b show the SC regime 
and c,d the regime just above T c . The upper and lower rows 
show scans along the a-axis (H, -1.5, -1.7) and 6-axis (1.5, 
K, 1.7), respectively. In order to obtain a meaningful colour 
representation, the intensity at 250 K was subtracted for E < 
38 meV and the data was corrected for a Q-linear background 
at all energies. At each individual energy, the colour scale 
was normalized to the peak intensity of the scan, allowing 
a better comparison of the Q-extent at different energies. 
The final wave-vector was fixed to 2.66 A below 38 meV 
and to 4.5 A above. Scans taken at the overlapping en- 
ergy 38 meV were used to bring both energy ranges to the 
same scale. Crossings of black lines represent measured data 
points. White lines connect the fitted peak positions of the 
constant-energy cuts. Dotted lines represent upper bounds 
on the incommensurability. 
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FIG. 3: Colour representation of the in-plane magnetic inten- 
sity collected at the TOF spectrometer at 5 K. The incident 
energy was fixed to Ei = 120 meV and the counting time was 
48 hours. In a intensity was integrated over the energy range 
E = (30 ± 3.5) meV and in b over E = (63 ± 5) meV. The 
data represent acoustic excitations and were collected close to 
the maximum L-component at L — 2.1 and 4.7, respectively 
(see Methods). A background quadratic in Q was subtracted 
and data are shown in arbitrary units. 



tary Notes and Ref. [20(). In contrast, at energies above 
E* the spectral weight is only moderately reduced upon 
heating above T* (Fig. la,b). We can thus distinguish 
a third temperature regime above T* , with a magnetic 
excitation spectrum differing distinctly from the spectra 
in the SC state and just above T c . 

Two aspects of our data in the PG state merit fur- 
ther comment. First, we observe the low-energy mag- 
netic spectral weight to increase upon cooling below T*, 
in contrast to the loss of spectral weight indicated by 
nuclear magnetic resonance (NMR) data [H. However, 
one needs to keep in mind that the energy scales probed 
by NMR and neutron scattering differ by three order of 
magnitude. The behaviour we observe can be viewed as 
a pile-up of spectral weight above a spin-pseudogap of 
« 15 meV, which is much larger than the NMR energy 
scale but consistent with prior neutron scattering work 
on underdoped YBa2Cu306+£c |2l|. Second, the magni- 
tude of the spin-PG implied by our data is much lower 
than the charge-PG observed in infrared spectroscopy [l[. 
The divergence of gap features in spin and charge chan- 
nels is not unexpected in a correlated metallic state close 
to the Mott insulator (where the spin gap vanishes, and 
the charge gap is ~ 2 eV). 

We now discuss the implications of our observations 
for the microscopic description of the PG state. Many 
features of the magnetic spectrum in the SC state are 
well described by models that regard the downward- 
dispersing branch as an excitonic mode in the spin-triplet 
channel below the SC energy gap 14|, 12J, |25j. If the 
downward dispersion is a manifestation of the d-wave 
symmetry of the gap, it is expected to disappear in 
the normal state, as observed. However, our observa- 
tion of coherent, highly anisotropic spin excitations for 
T c < T < T* requires a fresh theoretical approach. 
Specifically, the striking difference between the magnetic 



dynamics in the SC and PG states implies that these 
excitations cannot simply be regarded as incoherent pre- 
cursors of the resonant excitations below T c . Rather, 
they appear to be a signature of a many-body state that 
competes with superconductivity. 

The strong in-plane anisotropy of the spin dynamics in 
the PG state offers intriguing clues to the nature of this 
state. The fourfold orientational symmetry of the Cu02 
planes is broken in YBa2Cu3 0g+x by the CuO chains 
running along the 6-direction. However, weak-coupling 
calculations assuming a Fermi liquid state above T c in 
conjunction with recent ARPES-measurements [2rj| indi- 
cate that the hybridization between electronic states on 
CuO chains and Cu02 layers is not large enough to ex- 
plain the in-plane anisotropy at the low-energy end of 
the magnetic spectrum in the SC state Based on 

these model calculations it is quite unlikely that 



band structure effects alone can explain the much larger 
in-plane anisotropy and the strong temperature depen- 
dence we have observed above T c . The low-energy spin 
dynamics in the PG state thus appears to reflect a sus- 
ceptibility of the two-dimensional electron system to a 
weak in-plane symmetry-breaking field (induced by the 
CuO chains) that is markedly larger than that of a weakly 
renormalized Fermi liquid. Two possibilities have been 
discussed recently: (i) a "striped" state with quasi-one- 
dimensional order of spins and char ges that breaks ro- 
tational and translational symmetry 27|, l28| ; and (ii) a 
"nematic" or "Pomeranchuk" state that breaks the rota- 
tional symmetry only [B|, @, 0, S| ■ Experimental evidence 
for rotational symmetry breaking has also been gleaned 
from transport experiments 2!) . 

A static striped state can be ruled out for 
YBa 2 Cu306.6 based on the two-dimensional geometry of 
the high-energy spin excitations (Fig. 3), and on the ab- 
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FIG. 4: Temperature dependence of the peak intensity at 
three different positions in energy-momentum space sketched 
in reference to the "hour glass" dispersion in the SC state. 
In a, the energy transfer was fixed to E = 37 meV and the 
wave-vector transfer to Q = (0.5, 1.5, 1.7). In b, E — 30 meV 
and Q = (1.5,0.41,1.7) and in c, E = 30 meV and Q = 
(0.5, 1.5, 1.7). A weighted average of the background mea- 
sured at two representative points was subtracted. The units 
in the three panels were scaled to approximately the same 
intensity at T c . Error-bars indicate the statistical error. 
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sence of Bragg reflections indicative of static spin and 
charge order. However, predictions for the spin dynam- 
ics of fluctuating stripe 17 1 and Pomeranchuk states [l9j 
bear some resemblance to our data. Specifically, fluctuat- 
ing stripe segments with correlations lengths of the order 
of a few lattice spacings can result in an approximately 
"vertical", highly anisotropic dispersion surface at low 
energies, while the high-energy excitations are isotropic, 
as observed. However, calculations that explicitly con- 
sider the influence of superconductivity on the spin exci- 
tations of striped states [30] show only a rather modest 
effect (namely, the opening of a spin gap), in disagree- 
ment with our observations. Calculations for itinerant 
electron systems close to a Pomeranchuk instability [l|| 
are also in good overall agreement with our data; however 
in the SC state, the low-energy anisotropy of the spec- 
tral weight is significantly weaker than observed. Despite 
these promising first steps, an explanation of the unusual 
spin dynamics we have observed therefore remains an in- 
teresting challenge for theory. Its relationship to other 
fingerprints of the PG state, such as the presence of 
vortex-like excitations [9j and unconventional magnetic 
order [l(J, is also not understood at present. 

The full in-plane momentum resolution we have now 
achieved also opens a new window on the universality 
of the spin excitations in the cuprates. The high-energy 
spin excitations with four peaks along the diagonals of 
the Cu02 plaquettes are indeed common to all layered 
cuprate compounds thus far investigated. These exci- 
tations do not show marked changes at T e , and their 
intensity distribution exhibits at most a weak in-plane 
anisotropy. Our data show that the geometry and spec- 
tral weight distribution at lower energies depends more 
strongly on materials and temperature. In particular, 
we have shown that the spectra in the normal and su- 
perconducting states are qualitatively different. These 
differences may reflect the competition between high- 
temperature superconductivity and other ordering phe- 
nomena such as "nematic" or "Pomeranchuk" states. 

Correspondence and requests for materials should be 
addressed to B.K. (b.keimer@fkf.mpg.de). 
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Methods 

The experiments were performed on an array of 180 
individually detwinned YBa2Cu30g.6 crystal with super- 
conducting transition temperatures (midpoint) of T c sa 
61 K and width AT r » 2 K, determined for each crystal 
by magnetometry [13| . They were co-aligned on three 
Al-plates with a mosaicity of < 1.2°. The volume of the 
entire array was ~ 450 mm 3 , and the twin domain pop- 
ulation ratio was 94:6. 

Triple-axis measurements (Figs. 1, 2 and 3) were per- 
formed at the IN8 spectrometer at the Institut Laue 
Langevin (Grenoble, France) and the 2T spectrometer at 
the Laboratoire Leon Brillouin (Saclay, France). At both 
instruments, horizontally and vertically focusing crystals 
of pyrolytic graphite, set for the (002) reflection, were 
used to monochromate and analyze the neutron beam. 
Scans along a* and b* were carried out under identical 
instrumental resolution conditions by working in two dif- 
ferent Brillouin zones, (H, -1.5, -1.7) and (1.5, K, 1.7). 
No collimators were used in order to maximise the neu- 
tron flux, and graphite filters extinguished higher order 
contamination of the neutron beam. Extensive resolu- 
tion calculations confirm that the effects discussed in the 
main text are no resolution effects. Specifically, the spec- 
tra at 70K exhibit only structures that are significantly 
broader than the resolution function, so that the observed 
spectrum is almost undistorted by resolution effects. 

Time-of-flight measurements (Fig. 3) were performed 
at the MAPS-spectrometer of the ISIS spallation source, 
Rutherford Appleton Laboratory (Chilton, UK). The 
source proton current was 170 /iA and the Fermi chopper 
rotation frequency was set to 250 Hz. 

The wave vector is quoted in units of the reciprocal 
lattice vectors a*, b* and c* where a — 2ir/a* — 3.82 A, 
b = 3.87 A and c = 11.7 A. We choose its out-of-plane 
component Lq = 1.7 x (2n+ 1), n integer, to probe mag- 
netic excitations that are odd under the exchange of two 
layers within a bilayer unit. As even excitations show a 
gap of ps 55 meV and are much less T-dependent, they 
are presented elsewhere (25| . The incommensurability 
5 is defined as the deviation of the peak position from 
Qaf = (0.5,0.5). Even in the case, where two incom- 
mensurate peaks are close together or broad ("flat-top" 
structure), S can be obtained by fitting Gaussians to the 
data. 
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Supplementary material 




Figure SI: Colour representation of the magnetic intensity. 
Panels a,b show the SC regime, c,d the regime just above 
T c , and e,f the difference between both spectra. The upper 
and lower rows show scans along the a*-axis (H , -1.5, -1.7) 
and 6*-axis (1.5, K, 1.7), respectively. In order to obtain a 

meaningful colour representation, the intensity at 250 K was 

subtracted for E < 38 meV and the data was corrected for a 
Q-linear background at all energies. The final wave vector 
was fixed to 2.66 A below 38 meV and to 4.5 A 1 above. 
In contrast to Fig. 2 of the manuscript, where the colour 
scale was normalized to the peak intensity of the scan at 

each individual energy, we use an energy-independent colour 
scale here. Scans taken at the overlapping energy 38 meV 
were used to normalize both data sets. Crossings of black 
lines represent measured data points. White lines connect 

the fitted peak positions of the constant-energy cuts. Dotted 
lines represent upper bounds on the incommensurability. 

In Fig. SI, the data in the superconducting (5 K) 
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and pseudogap (70 K) states (same as in Fig. 2 of the 
manuscript) are reproduced in a different presentation 
scheme with a global colour normalization. The differ- 
ence spectrum 7(5 K) — 7(70 K) is also shown. The dif- 
ference signal comprises the downward-dispersing branch 
below the resonance energy E Tes = 37.5 meV, which 
draws its spectral weight from a limited range above and 
a more extended range below E Yes (negative signal in 
Fig. SI e,f). It is significantly less anisotropic than the 
PG spectrum itself (Fig. SI c,d), with respect to both 
the incommensurability S and the spectral weight distri- 
bution. Notably, the difference spectrum is very similar 
to its analog in almost optimally doped YBa2Cu3 06.85, 
which was shown to exhibit a nearly circular geometry in 
V. Hinkov et al., Nature 430, 650. 

Fig. S2 highlights again the difference between the 
"hour glass" and "vertical" dispersion relations in the 
SC and PG states, respectively. In the SC state, the pro- 
file along a* exhibits a single peak at E les = 37.5 meV, 
corresponding to the neck of the "hour glass" . In the PG 
state, the profile at this energy consists of two incom- 
mensurate peaks. The opposite trend holds along b* at 
energies below E les . The profile along b* is very narrow 
in the PG state, and incommensurate reflections cannot 
be resolved (dashed lines in Fig. 3 of the manuscript). 
Below E res , the incommensurability increases in the SC 
state, where the anisotropy between the dispersions along 
a* and b* is significantly reduced. 

Fig. S3 shows low-energy profiles deep in the SC state, 
along with data immediately below and above the SC 
transition temperature. While the scans at 5 and 50 
K show clear incommensurate peaks are identical within 
the experimental error, the profile in the scan at 70 K 
is flat-top, and the peak intensity is considerably re- 
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Figure S2: Comparison of constant-energy profiles at 37 
meV along a* , panel a, and at 33.5 meV along b* , panel b, 
in the SC and PG states. The data was taken from Figs, le 
and h of the manuscript. At 37 meV, the amplitudes of 
both scans were scaled to the same amplitude in order to 
compare their Q-width. Panels c and d show corresponding 
sketches of the dispersion relations in the SC state (grey) 
and the PG state (red). 
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Figure S3: Constant-i? scans at 31 meV along a* at different 
temperatures. For better comparison, a background 
measured at room temperature has been subtracted. 



duced. The data are complementary to the temperature 
sweeps at constant wave vector presented in Fig. 4 of the 
manuscript and demonstrate the abrupt modification of 
the magnetic spectrum at the SC transition temperature. 

Fig. S4 provides a synopsis of the spectra in the SC 
state at 5 K, the PG state at 70 K, and the normal state 
at room temperature, integrated along the wave vector 
component along a* . Whereas all three spectra are simi- 
lar at the highest energies, they exhibit qualitative differ- 
ences at lower energies, as described in the manuscript. 
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Figure S4: Energy dependence of the integrated intensity 
along a* in the SC state (5 K), the PG state (70 K) and the 
normal state (290 K). 



